Phases and phase transitions of hard-core bosons on a triangular ladder 
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We study hardcore bosons on a triangular ladder at half filling in the presence of a frustrating 
hopping term and a competing nearest neighbor repulsion V which promotes crystallization. Using 
the finite-size density-matrix renormalization group method, we obtain the phase diagram which 
contains three phases: a uniform superfluid (SF), an insulating charge density wave (CDW) crystal 
and a bond ordered insulator (BO). We find that the transitions from SF to CDW and SF to BO 
are continuous in nature, with critical exponents varying continously along the phase boundaries, 
while the transition from CDW to BO is found to be first order. The phase diagram is also found to 
contain an exactly solvable Majumdar Ghosh point, and re-entrant SF to CDW phase transitions. 

PACS numbers: 75.40.Gb, 67.85.-d, 71.27. +a 



I. INTRODUCTION 



Quantum phase transitions in low dimensional systems 
are of great interest because quantum fluctuations play a 
greater role in their physics than in their higher dimen- 
sional counterparts For instance in one dimensional 
systems, quantum fluctuations can inhibit the appear- 
ance of long range order that is seen in higher dimen- 
sional systems with the same symmetries. An example 
of this is the one dimensional spin 1/2 Heisenberg chain 
with nearest-neighbor antifcrromagnctic coupling which 
does not display any long range Neel order unlike its two 
dimensional square lattice counterpart. While quantum 
fluctuations in high dimensions maybe weak, geometrical 
frustration could prevent the occurrence of long range or- 
der. A naive expectation might be that quantum fluctua- 
tions and geometrical frustration, when present together 
would reinforce each other and severely inhibit ordering. 
However, there are several examples where the combina- 
tion of frustration and quantum fluctuations induces in- 
teresting types of order such as through order-by-disordcr 
transitions^ In one dimension, where quantum fluctua- 
tions are especially strong, their interplay with frustra- 
tion can give rise to interesting phases and phase transi- 
tions such as the recently proposed chiral Mott insulator 
in frustrated Bose-Hubbard ladders, a phase with a gap 
to all excitations and a staggered pattern of equilibrium 
currents 4^ 

In this paper, we study a frustrated one-dimensional 
system of hardcore bosons. In this model a boson on 
a site can hop to a neighboring site (with amplitude t) 
and also feels a density-density interaction (of strength 
V) from occupied neighboring sites. In addition to these 
terms, there is also a frustrating next nearest neighbor 
hopping term (of amplitude t'). The Hamiltonian that 



describes this t-t'-V model is 

H = -^(ata i+1 +h.c.)-*'^(ata i+2 +h.c.) 



n-i+i 



(1) 



where a] and a i are creation and annihilation operators 
for hard core bosons at site i, and rii = a\a i is the boson 
number operator at site i. Here we have the constraint 
that a\ 2 = a 2 =0, which avoids multiple occupancies 
of the lattice sites. Frustration in this model arises from 
taking t > and t' < 0. Further V > so that the in- 
teraction is repulsive. The model can be thought of as a 
zig-zag ladder of the sort shown in Fig.[T]with the nearest- 
neighbor and next-nearest neighbor hopping arising from 
the motion of the bosons between and along the legs of 
the ladder respectively. Note that this model does not 
have a simple representation in terms of spinless fermions 
due to the presence of the next-nearest-neighbour hop- 
ping term, which when expressed in terms of spinless 
fermions will not correspond to a simple hopping term. 
Thus, our model even with V = is non-trivial and does 
not have a ground state corresponding to a filled Fermi 
sea. In fact, as we will sec, for V = 0, aside from the 
trivial point t' — 0, there is only one other point corre- 
sponding to t' = —t/2, where the ground state can be 
obtained exactly. 

Eqn. [T] can be mapped onto a spin 1/2 Hamiltonian 
by identifying = a], = a 2 -, and Sf = (n, — 1/2). 
Under this transformation the Hamiltonian (JXJ) takes the 
form 



2t(s:s^ 1 + sfsf +1 ) + vs!s: +1 



-2t'(SfSf +2 + SfSf +2 ) 



(2) 
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FIG. 1: (Color online) Zig-zag ladder system representing 
the t — t' — V model in one dimension. The arrows are the 
representation of the hopping directions, t > and if < 0. 



By making a spin rotation on one sublattice, to change 
the sign of the next neighbor exchange coupling, Eqn. [5] 
can be recast as an anisotropic (XXZ type) spin-i model 
with a next-nearest-neighbor coupling. This model may 
then be viewed as a specific easy-plane deformation of 
the well-known SU(2) symmetric J\-Ji model which 
has been studied extensively^ - — . The SU(2) symmet- 
ric J\-Ji model exhibits gapless to gapped transition at 
J 2 « 0.241167^ Interestingly, when J 2 = Ji/2, the 
quantum system dimerizes in the presence of frustration 
and is described by the Majumdar- Ghosh (MG) modelii. 
Variants of the J\-J 2 model with spin anisotropy have 
also been studied. The earliest such study was by Hal- 
dane who added an anisotropy parameter to the nearest- 
neighbor coupling term while the next-nearest-neighbor 
term was left isotropic— More recently, a variant of 
the J\-Ji model with the same anisotropy in both the 
nearest-neighbor and next-nearest-neighbor terms has 
been studiedi 13 i 14 Our model given by Eqn. [2j while 
anisotropic in spin is different from both the anisotropic 
models mentioned above in that, it has only a single 
SU(2) symmetric point in parameter space (correspond- 
ing to V = 2t and if = 0) while the other two models have 
lines of points with SU(2) symmetry. Some of the salient 
features of our model like continuously varying critical 
exponents along the phase boundary and re-entrant tran- 
sitions are presumably due to the absence of such SU(2) 
symmetric lines, which will be discussed later. We will 
show that our model also has a point in parameter space 
analogous to the MG model when V = and t' = t/2, 
but without SU (2) symmetry, where we can obtain the 
ground state exactly. The model at this point can be 
thought of as a U(l) analog of the MG model. 

In this paper, we study the model given by Eqns. [T]at 
half filling numerically using the Density Matrix Renor- 
malization Group (DMRG) algorithm in the entire pa- 
rameter space of if < and V > with t > 0. We 
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FIG. 2: (Color online) Phase diagram of the half-filled t-t'-V 
model. There exist two gapless to gapped transitions, from 
SF to BO and SF to CDW for small \t'\ and V values, which 
are of the BKT type. There is also a direct transition from 
gapped BO to gapped CDW which is first order for large 
values of the interaction. The point at |t'| = 0.5 and V = 
is the Majumdar-Ghosh point and the dashed line is the 
theoretically calculated instability line for the CDW state. 
The hopping amplitude is set to t = 1. 

find that there are three phases, a superfluid (SF), bond- 
ordered (BO) phase and a charge density-wave (CDW). 
These are analogs of the gapless spin liquid, gapped 
dimerized phase and gapped Neel phase of the analo- 
gous anisotropic spin models. Our study mainly focusses 
on the phase boundaries and we show that the phase 
transitions in this model are of the Berezinski-Kosterlitz- 
Thoulcss (BKT) type and first order. Our results are 
shown in the phase diagram of Fig. [21 We find that 
along the lines of BKt transitions, there is a continuously 
varying critical exponent (corresponding to the Luttingcr 
paramter). Another interesting feature of the phase di- 
agram is a re-entrant SF-CDW-SF transition that exists 
in a part of parameter space and whose origin can be 
understood in terms of an instability of the CDW state. 

The remainder of the paper is organized as follows. 
In section II, we study the MG point analytically and 
obtain the exact ground state to better understand the 
BO phase and in section III, we study the CDW phase 
and obtain an analytical form for the phase boundary 
between this phase and the SF highlighting the reentrant 
phase transition. In section IV, we describe the details 
of our DMRG study and give expressions for the various 
quantities that have been calculated and in section V, we 
present our numerical results for the full phase diagram 
along with a discussion of its different features. 



II. MAJUMDAR-GHOSH POINT 

Consider the model H = — a\a i+1 +t' a]a i+2 +h.c.) 
for hard core bosons, where we have set the nearest neigh- 
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bor repulsion V = 0. We can rewrite this Hamiltonian in 
the language of spin-1/2 operators, which for t! = —t/2 
is given by 



while on the 4- sites of the CDW, we have 



MG 



Is 



-3/2+ (s%-sz +1 + s% +2 y 



Let us focus on a local state defined on 3 successive sites 
UJ+1J+2), of the form \ip) = \(tjl j+ i + ijtj+l^i+i), 
where (Jj+2 is an arbitrary spin state at site j + 2. It is 
easy to show that 

(S? - SJ +1 + SJ +2 )\i>) = ^Kt^+i + ^t 3 -+i)^ + a> (3) 

where (Tj+2 denotes the original spin at site j + 1 be- 
ing flipped. This means (5? - +S'^ f2 ) 2 |^} = dis- 
similarly, 

(5j'-5jV 1 +5jV 2 )|^) = ! ^l(t J l, + i+l J t, + i)^ +2 ), (4) 

so that (5V -5V + 1 +5j' +2 ) 2 |V) = ||# Thus, the product 
state 



Itili+i + liti+i) 



(5) 



j£even 



is an eigenstate of the Hamiltonian H^f Gl with energy 
per site |(— 3/2 + 1/2) = —t/2. Since a sum of three 



spin-1/2 operators, such as S y , 



qv 



i+Sj +2 i cannot have 



a minimum eigenvalue of magnitude less than 1/2, the 
state we have found is evidently a ground state. Another 
ground state, which is related by symmetry to the above 
ground state is simply 



II Itr^+i+VN+i)- 

j£odd 



(6) 



While we have not proved that these are the only two 
ground states of the Hamiltonian, our DMRG numerics 
indicate that this is the case. Note that the ground states 
obtained above are not products of singlets like the usual 
MG state and are thus not spin rotation invariant. This 
is not surprising since our Hamiltonian is also not spin 
rotation invariant unlike the J\-Ji model. 



III. PHASE BOUNDARY OF THE CDW STATE 
AND THE SF-CDW-SF RE-ENTRANT PHASE 
TRANSITION 

Consider the CDW state, which in spin language may 
be denoted as |... titi •■■)• Let us define Holstcin- 
Primakoff bosons h, such that on the t-sites of the CDW, 
we have^ 



1 



)\Sf 



,t 



(7) 



St = {h\h i -\);St = h\;Sr =h, 



(8) 



To quadratic order in the /i-bosons, we find the Hamil- 
tonian takes the form 

H = -t^hthl. + h.h^ + V^A 

i i 

- t'J^ihlh^ + hl^) (9) 

i 

Going to momentum space, and defining ^\ = (h,, h_ k ), 
we find the Hamiltonian 

TT xt (V — 2t' cos 2k —2tcosk \ T 

ff = E*U -2tcosfc V-2fcoB2jfe **' ( 10 ) 



k>0 



with eigenenergy Afe = y/ (V — 2t' cos 2k) 2 — (2t cos k) 2 . 
For large V, the spectrum has a gap ~ V. As V de- 
creases, however, the spectrum develops a soft mode sig- 
nalling an instability of the CDW state. For \t'\ < t/i, 
the instability develops at k = and k = tt below a crit- 
ical repulsive interaction V Ct i = 2(t — |t'|). For \t'\ > </4, 
the instability develops at an incommensurate wavevec- 
tors, fco and n — fco, where fco = cos _1 ( jrtt), below a 

critical repulsive interaction V c ^ = 2\t'\ + -jypy. The in- 
stability line is thus non-monotonic, with a minimum 
at \t'\ = and it approximately follows the phase 

boundary of the CDW state found numerically using the 
DMRG. The non-monotonicity of the instability line is 
responsible for the SF-CDW-SF re-entrant phase transi- 
tion as can be seen from Fig. [5] 



IV. DMRG TREATMENT OF THE MODEL 

We study the ground state of the model described by 
Eqn. [T] using the finite-size DMRG method with open 
boundary conditions! 16 ' 17 This method is best suited for 
(quasi-)one-dimensional problems*^ For most of our cal- 
culations we study system sizes up to 300 sites and retain 
128 density matrix eigen states with weight of the dis- 
carded states in the density matrix less than 10 -6 . When 
t! = 0, the model can be solved exactly using the Bethe 
ansatai^ and there exists a transition from a gapless SF 
to a gapped CDW phase like (. . . 10 10 1010 . . .) phase 
at V = 2t. Here 1 (0) represents the presence(absence) of 
a boson at a particular site. However, as we will shown 
when the value of t' is finite, the phase diagram is much 
richer. In the present work we calculate various different 
physical quantities to characterize the phases and phase 
transitions of our model. The resulting phase diagram is 
shown in Fig. [2] 

We now list the various quantities we have calculated 
to identify the phases in our model. These quanti- 
ties have been used previously to identify similar phases 
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in related models*^ In order to distinguish between 
the gapped and gapless phases we calculate the single- 
particle excitation gap 

G L = E(L, N + l) + E(L, N - 1) - 2E(L, N). (11) 

In Eq. (|TT1) . E(L, N) is the ground-state energy of a sys- 
tem with L sites and N bosons. 

The CDW phase and the transition into it can be 
studied by calculating the structure factor, which is the 
Fourier transform of the density-density correlation func- 
tion 

= £2 £ eifc(i ~ J) (<^> - (12) 

The BO phase has a non-zero value of the bond-order 
parameter 

i 

where 

Bi={c\c i+1 +c\ +rCi ), (14) 

and this is thus the quantity we calculate to identify this- 
phase and the phase transition into it. 

The commensurate to incommensurate transition can 
be tracked by identifying the k vector for ordering in a 
phase. This vector can be found by looking for a peak in 
the momentum distribution function 

Finally, to characterize the phase boundaries, we calcu- 
late the correlation function 

r(r) - (ajar), (16) 

which as we will see can be used to define the Luttinger 
parameter along the phase boundaries. In the remainder 
of the paper, we set t = 1. 

V. RESULTS AND DISCUSSION 

In the first part of this section we discuss the results of 
our calculation for V — 0. When t' = 0, the system is a 
gapless SF with finite momentum distribution at k = 0. 
As \t'\ increases, the system enters a gapped phase. The 
gapless to gapped transition can be seen by calculating 
the single particle excitation gap as given in Eqn. 1111 
We plot the extrapolated gap Gl->-oo as a function of \t'\ 
in Fig( It is clear that there exists a transition to 
a gapped phase at \t'\ ~ 0.35. The gap arises because 
of frustration which tries to make the ground state have 
bond order (BO) of the sort described in our analysis 
of the MG point. The bond energy is different for odd 




FIG. 3: (Color online) The extrapolated gap Gl — > oo plotted 
as a function of \t'\ for V = 0. The gapless to gapped phase 
transition occurs at \t'\ m 0.35 



and even bonds in this phase. This emergence of the 
BO phase can be tracked by calculating the BO order 
parameter (Obo) as given in Eqn. 1131 In the BO phase, 
Obo is finite in the thermodynamic limit and vanishes 
for small values of \t'\. We plot Obo as a function of 
\t'\ in Fig. H for three different lengths, L = 100,200 and 
300. The value of Obo is maximum at \t'\ « 5.9 and 
decreases on either side of the maximum. To the left of 
the maximum, i.e. when \t'\ decreases, the jump becomes 
sharper and sharper as the length increases and appears 
to be saturating to a small value as \t'\ decreases. At 
\t'\ = 0.5 the value of Obo = 0.5 and independent of 
length as expected for the MG point. 

Our analysis of sect. IIIII shows that there is a commen- 
surate to incommensurate transition at \t'\ = 0.25. We 
can verify this by calculating the momentum distribution 
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FIG. 4: (Color online) Obo is plotted as a function of \t'\ for 
L — 100, 200 and 300. The increases of the value of Obo as 
\t'\ increases indicates the transition to the BO phase. The 
vertical dashed line indicates the location of the MG point 
and the horizontal dashed line, the value of Obo at this point, 
which is 0.5 indpendent of system size. 
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n(k) as a function of k as shown in Fig. [5] Here we see 
that the momentum distribution has one maximum at 
k = up to about \t'\ « 0.5 after which the maxima shift 
to values of k ^ 0. This indicates the commensurate to 
incommensurate transition. The peak position shifts to- 
wards k = ±7r/2 as \t'\ — > oo where the two- leg ladder 
can be considered as a single chain. The fact that the 
transition does not occur at \t'\ = 0.25 can be attributed 
to the hard core nature of the bosons. The analysis of 
sect. IIIII assumed non-interacting particles. 

Now we turn our attention to V ^ 0. When \t'\ = it 
is known that at V = 2, the model exhibits a transition 
from SF to CDW phased. As the value of \t'\ is 
increased we get three different scenarios, 

(1) SF - CDW transition as a function of V for small 
values of \t'\ 

(2) BO - SF - CDW as a function of V for intermediate 
values of V. 

(3) BO - CDW transition as a function of V for large 
values of \t'\ 

It is to be noted that the SF phase is gapless where as 
the BO and CDW phases are gapped. When the value 
of \t'\ is small, the SF-CDW transition occurs at V < 2 
since the presence of negative \t'\ suppresses the effect 
of t. This allows the bosons to stabilize in the CDW 
ground state even for a small nearest neighbor interac- 
tion. This can be seen from Fig. [5] where the upper phase 
boundary represents the SF - CDW transition for small 
values of \t'\. As \t'\ increases the value of V at which the 
transition takes place decreases gradually till it reaches 
a minimum and then increases again as expected from 
the analysis of sect. IIIII This indicates a re-entrant phase 
transition where at fixed V , an increase of \t'\ will drive 



the system from SF to CDW and back to SF again. It 
is interesting to note that the lowest point of the phase 
boundary obtained from the numerics is not too far from 
that predicted from the analytical calculation of sect. IIIII 

The phase boundary between the SF and BO phases 
is also shown in Fig. [2j This boundary originates at V = 
and |t'| = 0.33 and moves upwards as shown. The 
supcrfluid region is pinched off by the approach of the 
SF-BO phase boundary towards the SF-CDW boundary. 
Crudely speaking, V prefers the formation of CDW order, 
\t'\ favors BO while the nearest neighbor hopping causes 
SF order. As long as both V and \t'\ are smaller than 
or close to 1, the nearest neighbor hopping term ensures 
that the transition from BO to CDW has an intervening 
region of SF. However, once V and \t'\ start becoming 
larger than 1, there is a direct transition from CDW to 
BO, which appears to be first order from our calculations. 

We now discuss our characterization of these transi- 
tions. It is known that for t' = 0, the transition from 
SF-CDW exhibits a BKT type scaling of the gap. The 
correlation function T(r) ~ l/r n with 77 = 1 at the criti- 
cal point with an antiferromagnetic modulation and also 
a log correction. We find a similar BKT type scaling of 
the gap along both the SF-CDW phase boundaries and 
SF-BO phase boundaries. The scaling of the gap at the 
BKT transition can be used to locate the critical point of 
the transition fairly accurately.—. The critical exponent 
77 can also be measured along the phase boundaries and 
we find this exponent to be varying continuously along 
both phase boundaries. We do not have the numerical 
accuracy to track this variation all the way up to the 
point where the boundaries appear to merge. It appears 
that there is a first phase transition from BO to CDW for 
larger values of \t'\ and V. In the following subsections, 
we present data supporting each of these claims. 





FIG. 5: (Color online) Momentum Distribution n(k) is plot- 
ted for different \t'\. There exists a commensurate to incom- 
mensurate transition at \t'\ > 0.5 which can be seen as the 
onset of two peaks at k 7^ 0, n. 



FIG. 6: (Color online) The scaled gap LG' L is plotted as a 
function of xl- (Inset) The scaled gap plotted as a function 
of \t'\ for V = 0. This shows the SF to BO transition at 
\t'\ = 0.33. 
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FIG. 7: (Color online) The scaled gap LG' L is plotted as a 
function of xl- (Inset) The scaled gap plotted as a function 
of \t'\. This shows the SF to BO transition at It'l = 0.32. 



FIG. 8: (Color online) The scaled gap LG' L is plotted as a 
function of xl- (Inset) The scaled gap plotted as a function 
of V. This shows the SF to CDW transition at V = 1.31. 



A. BKT scaling and gapless to gapped transition 

Since the SF phase is gapless and the BO and CDW 
phase are gapped, we can use scaling of the gap Gl given 
by Eqn. [IT] to locate the transition. Gl is computed for 
lattice up to 300 sites. At the BKT transition from SF 
to BO the gap closes as 



G 



exp 



j\\V\-\t'c\\ 



(17) 



where a is a constant. 

The correlation length £, which scales at the critical 
point as the inverse of the gap is finite in the gapped 
phase and diverges in the gapless SF phase. We use the 
following finitc-size-scaling relation for the gap in the re- 
gion close to the phase transition, 



LG L x 1 



1 



2 In L + C 



(18) 



where F is a scaling function and C is an unknown con- 
stant to be determined. In the region close to the criti- 
cal point and within the SF phase, the values of F(£/L) 
is expected to be system-size independent, i.e., plots of 
LG' L = LG L [1 + 1/ (2 InL + C)] as function of f for dif- 
ferent system sizes should intersect in that region. Also, 
the curves obtained by plotting LG' L as function of £/L 
for several values of L should be system-size independent. 
Therefore, the plots of LG' L as function of xl = In L — In £ 
for different lengths collapse in the critical region. We 
obtain the values of a, C, and \t' c \ for the best possible 
collapse of the data in the gapped side where the cor- 
relation length diverges as £ ~ exp a/J\\t'\ — \t' c \\ . A 

similar procedure can be used for the other part of the 
phase diagram i.e. the SF to CDW transition by replac- 
ing \t'\ by V. The accuracy of this method has been 



tested by locating the Heisenberg point at t' = 0J^. It 
is found to be at V = 2.02 ± 0.01, very close to the an- 
alytical result. Here, we also find the critical point for 
the SF to BO transition at \t'\ = 0.33 ± 0.01 for V = 
which is consistent with the value obtained previously in 
equivalent spin models.— 

We use the above technique to obtain the boundaries 
between gapped and gapless phases in the phase diagram 
Fig. [2j Another possible way to obtain the boundary is 
to extrapolate the gap to L — > oo and locate the points 
at which it goes from being non-zero to zero. However, 
we find that that the former technique is more accurate 
than the latter one and thus we use the BKT scaling form 
of the gap to locate the transition. In Fig. [7] we show the 
scaling of the gap along the SF - BO transition boundary. 
The collapse of the curves is obtained by plotting LG' L 
vs xl within the gapped phase for \t' c \ = 0.325±0.005 and 
V = 0.5 (main panel). The plots of LG' L vs \t'\ (insets), 
for V — 0.5 and for three values of L show that the 
LG' L curves intersect at the critical point \t'\ = 0.32. 
Similarly in Fig. [8J we show the collapse of the LG' L 
vs xl data for \t'\ = 0.3 when V = 1.31 (main panel) 
and the intersection of the LG' L vs V curves for different 
lengths at the critical point (inset). 

We now verify the locations of the SF-BO and SF- 
CDW phase boundaries by using the scaling of the BO 
order parameter and the density-density structure factor. 
In the BO phase Obo is finite and zero in the SF phase. 
The value of Obo is equal to 0.5 at the MG point. In 
order to see the the transition from BO to SF and then 
to CDW we start from the MG point i.e. from \t'\ = 0.5 
and then move along the V axis. In Fig. [5] we plot Obo 
as a fuction of V for L = 100, 200, 300. It can be seen 
that the value of Obo decreases as we increase V. The 
decrease is faster for large lengths decrease implying that 
Obo tends to zero in the thermodynamic limit. In or- 
der to see the actual transition point we perform a finite 
size scaling of Oso^r^ In Fig. [JJjl we plot Obo as a 



V 



FIG. 11: (Color online)Finite size scaling of S(n) shows the 
FIG. 9: (Color online)Oso is plotted as a function of V at SF-CDW transition. 
\t'\ = 0.5 for different lengths. 

B. Critical exponent across the phase boundary 




•-•v= 


1.1 


■-■ v= 


1.2 


♦-♦v= 


1.3 


A-AV= 


1.4 




1.5 


<-< v= 


1.6 


v= 


1.7 




1.8 


-*v= 


1.9 




2.0 



0.05 



0.15 



1/L 



0.5 



FIG. 10: (Color online) Obo is plotted as a function of 
l/L -0 ' 5 for different V at \t'\ = 0.5 showing the BO-SF tran- 
sition. 



function of L~ 5 for different values of V and then ex- 
trapolate to L — > oo. It is evident that the curve for 
V = 1.6 extrapolates to zero and the curves for V > 1.6 
extrapolate to finite values showing the transition to BO 
phase for values of V ~ 1.5. This result is in accordance 
with the phase diagram obtained by scaling of the gap 
where the transition point for the BO to SF transition is 
at V = 1.53. In order to understand the SF-CDW tran- 
sition we perform a similar scaling of the density-density 
structure factor as given in Eqn. 1121 An extrapolation 
shows that the <S(7r) for V > 1.6 tends to a finite value in 
the thermodynamic limit. However, for values of V < 1.6 
the curves appear to extrapolate to zero. This is also in 
accordance with the phase diagram where the SF-CDW 
transition occurs at V — 1.62 for \t'\ = 0.5 



We now compute the critical exponents across the two 
BKT phase boundaries representing the SF-CDW and 
SF-BO transition. This exponent 77 is defined by the 
relation 



r(r) ~ l/r r > 



(19) 



where r = \i — j\, i and j are the lattice indices, rj can 
thus be obtained by means of a straight line fit to T(r) 
as a function of r on a log-log scale. To avoid bound- 
ary effects, which can corrupt such a fit, we discard data 
obtained from the edges of our numerical system and 
use only data obtained from the bulk. The exponent 
r) obtained this way will have a dependence on system 
size and the number of DMRG states kept in the cal- 
culation, which we have investigated. The plots of T(r) 
vs. r arc shown in Fig. [12] for different points along the 
phase boundaries for L = 500 and 128 states. About 
50 sites from the bulk were used for the fit. The error 
bars mentioned are for the linear fits to the data. The 
values of rj while always close to 1, seem to be varying 
along parts of the phase boundaries. r\ was found to 
decrease for all points with increasing L and increasing 
number of DMRG states, so it is reasonable to believe 
that rj < 1 along sections of the phase boundary. We 
find r\ > 1 for some parts of the phase boundary such 
as t' = -0.3, V = 1.3 and t' = -0.35, V = 1.0 and it is 
likely that the value of 77 for these points will eventually 
go to 1 or lower with increasing system size and number 
of DMRG states which will have to be confirmed by more 
extensive numerical calculations. Further, the error in 77 
introduced due to the error in the determination of the 
critical points is minuscule and less than that due to the 
linear fits. 

As can be seen from Fig. Q21 77 appears to be increas- 
ing as we move to larger values of \t'\ along the SF-CDW 
boundary. It is known that 77 = 1 exactly at rj' = 
but from our data it appears that rj drops as soon as a 
non-zero t' is introduced and once again rises towards 1. 
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FIG. 12: (Color online)r(r) plotted as a function of r for 
10 < r < 50. The left panel is for the critical points along the 
SF-CDW boundary and the right panel is along the SF-BO 
boundry of the phase diagram shown in Fig([2} ■ The symbols 
are the value of T(r) and the solid lines are the fitted function 
of the form a/r v . The errors are obtained from the linear fit 
to the data. It can be clearly seen that the value of r\ increases 
by increasing the values of V and \t'\ as we move along the 
SF-CDW and SF-BO respectively. 



Along the other phase boundary, r\ starts from a value 
close to but less than 1 and increases as we move along it 
to larger values of V. It is likely that the two boundaries 
merge exactly at the point where rj = 1, which is also 
where the line of first-order transitions begins (a multi- 
critical point). Another possibility is that the two phase 
boundaries merge before the line of first order points be- 
gins and there is an intermediate section along the phase 
boundary between CDW and BO, where rj varies con- 
tinuously, i.e. the transition is Gaussian in nature. At 
the moment, we do not have the numerical accuracy to 
resolve these two scenarios. 

It is interesting to note that we obtain a value of r/ 
different from 1 and changing continuously at least along 
parts of the boundaries. This is not the case for previ- 
ously studied anisotropic spin models. It has been argued 
that in those models, the value of rj is pinned to 1 along 
both the SF-CDW and SF-BO phase boundaries!^. A 
rcnormalization group (RG) analysis of the sine-Gordon 
theory for these models shows that this is due to the fact 
that the zero umklapp line intersects a line of SU (2) in- 
variant points in the phase space of the models!^. The 
point of intersection happens to be the transition point 
from spin fluid (SF in the language of hardcore bosons) to 
dimer order (BO in the language of hard core bosons) and 
has rj = 1. The RG flow along both phase boundaries is 
towards this point thus pinning the exponent along them 
to 1. Further, in this RG analysis, the transition from 
CDW to BO is always continuous with a continuously 
varying rj > 1, a Gaussian transition. 

Our model does not have an SU (2) symmetric point 
separating the SF and BO phases or a line of SU (2) sym- 
metric points in parameter space and thus will not have 



0.03 




FIG. 13: (Color online)5'(7r) is plotted as a function of V 
for \t'\ = 1.0 showing the BO-CDW transition. (Inset)The 
— ¥ oo is plotted as a function of V to confirm the first 
order nature of the BO-CDW transition. 

the same RG flow diagram. There is hence no reason a 
priori to expect the rj value to be pinned to 1 along the 
phase boundaries. However, since the transition out of 
the SF phase might still be expected to be governed by a 
sine- Gordon theory in which umklapp is not relevant, rj 
should be less than or equal to 1 along these boundaries. 
Thus, it can be expected that the values of ?] > 1 we 
seem to obtain for certain points along the phase bound- 
aries will settle down to 1 or lower as the system size and 
number of DMRG states are increased. 

An analytical understanding of the phase boundaries 
and critical exponents will require a detailed field theo- 
retical study of the underlying sine-Gordon type action, 
which will be presented in a separate paper. 



C. Gapped to gapped phase transition 

In this section we study the transition from BO to 
CDW phase at large \t'\ values. For \t'\ > 0.6 the SF 
phase shrinks very slowly and appears to finally disap- 
pear at \t'\ as 0.7. After this point the BO phase slowly 
undergoes a direct transition to the CDW phase. To 
study the phase transition we use the CDW structure 
factor as an order parameter. We plot the extrapolated 
values of S(tt) as a function of V in Fig. [T2]for \t'\ = 1.0. 
The sudden jump in the value of S(tt) at V « 2.9 implies 
a first order transition from the BO to the CDW phase. 
This transition is further verified by examining the corre- 
lation length of the system. At the first order transition 
the single particle excitation gap remains finite. As men- 
tioned earlier, the correlation length £ oc 1/G. There- 
fore at the first order transition l/£ should remain finite. 
We show the extrapolated values of 1/£l->oo in the inset 
of Fig. [13] It is evident from the figure that while ap- 
proaching the transition from the BO side, 1/£l^oo(Gl) 
decreases rapidly, reaches a minimum at the transition 
point V = 2.97 and then increases as the system enters 
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FIG. 14: (Color online). The derivative of the ground state 
energy with respect to V as a function of V. A discontinuity 
appears in the first derivative, showing that the ground state 
energy has a kink as a function of V, indicating a discontin- 
uous transition. The "latent heat" , proportional to the mag- 
nitude of the discontinuity seems to increase with increasing 
\t'\. 

the CDW phase. The gap remains finite at the minimum 
implying the absence of an SF phase. A further indica- 
tor of the first order nature of the transition comes from 
looking at the derivative of the ground state energy with 
respect to V, dE/dV shown in Fig.[l4]for different values 
of \t'\. It can be seen that there is the appearance at a 
discontinuity in the derivative, which seems to get more 
pronounced with \t'\. This shows that the transition is 
discontinuous and the "latent heat" associated with it 
increases as one moves up the phase boundary. 

The sine-Gordon theory for the CDW to BO transition 
predicts a continuous transition with 77 > 1 varying con- 
tinuously along the phase boundary^. However, it has 
been pointed out that this analysis ignores higher order 
umklapp terms, which can drive the transition first or- 
der, which is what we appear to be seeing herein. We 
emphasize, once again, that we do not have the resolu- 
tion to clearly say if the first order line begins from the 
point of intersection of the SF-CDW and SF-BO phase 



boundaries. There could be a small section of the phase 
boundary where the Gaussian transition between CDW 
and BO is seen. However, our numerics suggest that for 
sufficiently large values of V along the CDW-BO phase 
boundary, the transition is of first order. 



VI. CONCLUSION 

We have presented a detailed study of the hardcore 
boson in a one dimensional lattice in the presence of 
frustrated next-nearest neighbour hopping and the near- 
est neighbour interaction using the finite size DMRG 
method. The ground state phase diagram has three 
phases, SF, CDW and BO with continuous transitions 
along the SF-CDW and SF-BO phase boundaries and 
first order transitions along the CDW-BO phase bound- 
ary. The SF-CDW phase boundary is not monotonic giv- 
ing rise to a re-entrant phase transition. Further, the 
critical exponent 77 is different from 1 and varies contin- 
uously along the SF-CDW and SF-BO phase boundaries 
in contrast to other anisotropic frustrated models. 
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